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Abstract 

We calculate the 6-dependent self-energy of the top quark at 0(g 2 a s ) by 
using a general massive two-loop algorithm proposed in a previous article. From 
this we derive by unitarity the 0(a s ) radiative corrections to the decay width 
of the top quark, where all effects associated with the b quark mass are included 
without resorting to a mass expansion. Our results agree with the analytical 
results available for the 0(a s ) correction to the top quark width. 
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Precision measurements of the electroweak parameters are a powerful tool for test- 
ing the validity of the standard model and searching for new physics. The impressive 
accuracy attained experimentally, which is also expected to improve in the future, 
makes a complete two-loop analysis of the data necessary. 

However, on the theoretical side, such a complete two-loop electroweak analysis 
is far from being a simple task. Apart from the proliferation of diagrams, a special 
kind of technical problem is encountered in electroweak two-loop calculations. These 
processes involve particles with different masses, and in general need to be evaluated 
at finite external momenta. It has been known already for some time that general 
massive two-loop Feynman diagrams, when evaluated at nonvanishing external mo- 
menta, lead to complicated and often unknown special functions. See, for instance, 
ref. [|I| which relates the sunset two-point topology to the Lauricella function. 

Recent progress in two-loop electroweak analyses was mainly attained by using 
mass or momentum expansion methods. We note that in some situations mass expan- 
sions turn out to converge well, such as the top mass expansion of the a s corrections 
to the Z —>■ bb process |2[] , while in other cases subleading terms are known to be sub- 
stantial || . In certain situations it was possible to recover the exact function starting 
from an expansion Q. On the other hand, when the process under consideration 
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involves more than one ratio of masses, recovering the exact result from an expansion 
would be difficult, the functions involved being complicated. 

In order to circumvent the use of a mass or momentum expansion, we proposed in 
ref. || a general framework for treating two-loop Feynman graphs by a combination of 
analytical and numerical methods. Given a Feynman diagram, its tensorial structure 
is first reduced analytically into a combination of ten special functions hi, which 
are defined by fairly simple integral representations. The result of this analytical 
reduction is further integrated numerically. 

A subset of the ten special functions hi ||, namely hi and hi, are sufficient for 
treating two-loop corrections in a theory involving only scalars and no derivative 
couplings. This is the case of radiative corrections of enhanced electroweak strength 
in the standard model, and in ref. || this method was used for deriving corrections 
of leading power in the Higgs mass. The resulting radiative corrections agree with 
independent calculations which use different techniques 0. 

Going from a scalar theory to a full renormalizable theory with fermions and 
derivative couplings, such as the standard model, increases a lot the complexity of 
the problem. The analytical reduction of the tensor structure of a graph into hi 
functions along the lines of ref. || needs to be handled by a symbolic manipulation 
program such as FORM or Schoonschip, and results in rather lengthy expressions. 
Once the reduction is performed, the numerical integration methods are the same as 
those used in ref. ||, with the only difference that more computing power is needed 
to handle the number of hi functions which need to be integrated. 

In ref. [§] we have shown how these methods can be used for calculating momen- 
tum derivatives of two-loop two-point functions around a finite, non-zero value of the 
external momentum. Such momentum derivatives are encountered for instance when 
evaluating wave function renormalization constants on-shell. 

In this letter we show how this method can be used for calculating two-point 
functions of physical interest. We note that, due to the general nature of the tensor 
reduction algorithm given in ref. ||, the inclusion of more than two external legs 
proceeds in the same way. The difference is that the resulting expressions are more 
complicated than those stemming from two-point functions, and more computing 
power is needed for performing higher-dimensional numerical integration. 

In figure 1 we show the two-loop Feynman graphs relevant for the 6-dependent 
self-energy of the top quark at order g 2 a s . 

As for the counterterm structure, we only show in figure 1 the counterterm di- 
agrams which are necessary for renormalizing the imaginary part of the self-energy. 
These are the on-shell one-loop QCD counterterms of the 6-quark mass bm^ and of 
the top wave function renormalization constant 5Z t : 
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Figure 1: The two-loop Feynman graphs which contribute to the b-mass dependent 
correction of 0(a s g 2 ) to the top self-energy. Only the counterterm diagrams are 
shown which are needed for subtracting the infinities of the imaginary part of the 
self-energy, which gives the 0(a s ) correction to the t — > W + b decay. 
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Here, /i is the 't Hooft mass scale, and m g is the gluon mass infrared regulator. 
N c = 4/3 is the color factor. 

The imaginary part of the two-loop self-energy on-shell is physical. Writing the 
top self-energy as E(p -7) = Ei(p • 7) + 75 • S 75 (p ■ 7), the top decay width is given by 
F t = 2 • ImTji{p ■ 7 = m t ). This can be compared to known analytical and numerical 
results for the 0(a s ) QCD corrections to the top decay width 0, and thus provides 
a nontrivial test of the two-loop algorithm. 

The algebraic reduction of the tensorial structure of the graphs was done by im- 
plementing the reduction formulae of ref. into a symbolic manipulation program. 
The resulting intermediary decomposition into hi functions is too lengthy to repro- 
duce here. Instead, we give the final results obtained after numerical integration over 
the analytical decomposition. 

Here we would like to discuss briefly the treatment of infrared divergence associ- 
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Figure 2: The real and the imaginary parts of Y>i(p ■ 7 = m t ) at two-loop (see text), 
given by the graphs of figure 1, as a function of the top mass and of the gluon infrared 
mass regulator. An overall factor a s N c GFMwm t / V2/ (2ir) 8 is understood. 

ated with the massless gluon. To regularize the infrared singularities, we introduce 
a mass regulator for the gluon. We note that in higher-order QCD calculations in- 
troducing a gluon mass would affect the Slavnov- Taylor identities. In such a case, a 
different approach is needed to treat the infrared singular diagrams. One possibility is 
to try to extract analytically the IR singular part of the diagram and to treat this by 
dimensional regularization, while the remaining IR finite part of the diagram can be 
calculated by numerical integration. Whether or not this can always be done in a sim- 
ple way is beyond the scope of this article. We also note that our approach is mainly 
designed to treat massive graphs, while other approaches are available for massless 
QCD calculations. However, this being an 0(g 2 a s ) correction, the use of a gluon mass 
regulator is legitimate in this case. The physical quantity ImEi(p-^ = m t ), which is 
associated with the top quark width, is free of infrared singularities, and we checked 
numerically that indeed the result becomes independent of the gluon regulator when 
the mass regulator is much smaller than the b mass. 

In figure 2 we give numerical results for the ultraviolet finite part of the quantity 
• 7 = m t ) derived from figure 1; the inclusion of the counterterm contributions 
shown in figure 1 makes only the imaginary part (which is physical) of the self-energy 
finite, both UV and IR. We plot the results for a range of the top mass, and assume 
G F = 1.16637 • lCT 5 GeV~ 2 , m w = 80.41 GeV, m b = 4.7 GeV, and a s (m t ) = .108. 

In table 1 we give numerical values for the 0(a s ) correction to the top decay 
rate t — > W + b obtained from the imaginary part of the self-energy. Therefore, the 
0(a s ) QCD correction is obtained as an inclusive quantity, integrated over the gluon 
spectrum. To check the size of the finite mass of the b quark, we ran our programs 
in the vanishing b mass limit. At tree level the finite b mass amounts to a 3-4 MeV 
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m t [GeV] 


160 


165 


170 


175 


180 


T* ree [GeV] 


1.127 


1.260 


1.402 


1.553 


1.712 


ST i-loop [GeV j 


-.092 


-.104 


-.117 


-.132 


-.149 



Table 1: The 0(a s ) correction to the top decay t — > W + b as obtained from the 
imaginary part of the two- loop top self-energy of figure 1, integrated numerically. We 
took G F = 1.16637-10- 5 GeV" 2 , m w = 80.41 GeV, m b = 4.7 GeV, and a s (m t ) = .108. 

correction to the width, and in the 0(a s ) correction the effect is negligible. These 
results agree with existing calculations of 0(a s ) corrections to t — > W + b [[|, which 
provides a good test for the two-loop tensor decomposition and numerical integration 
algorithm. 

To conclude, we have shown that the general methods proposed in ref. f| can 
be used to calculate physical radiative corrections. We treated at two-loop order 
the 6-mass dependent self-energy of the top quark at 0(g 2 a s ). From the imaginary 
part of the self-energy we extract the 0(a s ) corrections to the top decay process 
t — > b+W, and find agreement with existing results for this process. The calculation is 
performed while respecting the actual mass and momentum kinematics of the process, 
and without resorting to a mass or momentum expansion of the diagrams. 
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